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Abstract: - The stability of an interconnected power system is its ability to return to normal or stable operation 
after having been subjected to some form of disturbance. With interconnected systems continually growing in 
size and extending over vast geographical regions, it is becoming increasingly more difficult to maintain 
synchronism between various parts of the power system. Modern power systems have many interconnected 
generating stations, each with several generators and many loads. So our study is focused on multi -machine 
stability . 
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I. INTRODUCTION 

Successful operation of a power system depends largely on the engineer's ability to provide reliable and 
uninterrupted service to the loads. The reliability of the power supply implies much more than merely being 
available. Ideally, the loads must be fed at constant voltage and frequency at all times. The first requirement of 
reliable service is to keep the synchronous generators running in parallel and with adequate capacity to meet the 
load demand. Synchronous machines do not easily fall out of step under normal conditions. If a machine tends 
to speed up or slow down, synchronizing forces tend to keep it in step. Conditions do arise, however, such as a 
fault on the network, failure in a piece of equipment , sudden application of a major load such as a steel mill, or 
loss of a line or generating unit., in which operation is such that the synchronizing forces for one or more 
machines may not be adequate, and small impacts in the system may cause these machines to lose synchronism. 

A second requirement of reliable electrical service is to maintain the integrity of the power network. 
The high- voltage transmission system connects the generating stations and the load centers. Interruptions in this 
network may hinder the flow of power to the load. This usually requires a study of large geographical areas 
since almost all power systems are interconnected with neighboring systems. 

Random changes in load are taking place at all times, with subsequent adjustments of generation. 
One may look at any of these as a change from one equilibrium state to another. Synchronism frequently may be 
lost in that transition period, or growing oscillations may occur over a transmission line, eventually leading to its 
tripping. These problems must be studied by the power system engineer and fall under the heading "power 
system stability". 

II. STUDY OF SWING EQUATION 

2.1 STABILITY 

The tendency of a power system to develop restoring forces equal to or greater than the disturbing forces to 
maintain the state of equilibrium is known as " STABILITY'. 

The problem of interest is one where a power system operating under a steady load condition is 
perturbed, causing the readjustment of the voltage angles of the synchronous machines. If such an occurrence 
creates an unbalance between the system generation and load, it results in the establishment of a new steady- 
state operating condition, with the subsequent adjustment of the voltage angles. The perturbation could be a 
major disturbance such as the loss of a generator, a fault or the loss of a line, or a combination of such events. It 
could also be a small load or random load changes occurring under normal operating conditions. Adjustment to 
the new operating condition is called the transient period. The system behavior during this time is called the 
dynamic system performance, which is of concern in defining system stability. The main criterion for stability is 
that the synchronous machines maintain synchronism at the end of the transient period. 

So we can say that if the oscillatory response of a power system during the transient period following a 
disturbance is damped and the system settles in a finite time to a new steady operating condition, we say the 
system is stable. If the system is not stable, it is considered unstable. This primitive definition of stability 
requires that the system oscillations be damped. This condition is sometimes called asymptotic stability and 
means that the system contains inherent forces that tend to reduce oscillations. This is a desirable feature in 
many systems and is considered necessary for power systems. The definition also excludes continuous 
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oscillation from the family of stable systems, although oscillators are stable in a mathematical sense. The reason 
is practical since a continually oscillating system would be undesirable for both the supplier and the user of 
electric power. Hence the definition describes a practical specification for an acceptable operating condition. 
The stability problem is concerned with the behavior of the synchronous machines after a disturbance. For 
convenience of analysis, stability problems are generally divided into two major categories -steady state stability 
and transient state stability and transient state stability. 

2.2 SWING EQUATION 

Under normal operating conditions, the relative position of the rotor axis and the resultant magnetic 
field axis is fixed. The angle between the two is known as the power angle or torque angle. During any 
disturbance, rotor will decelerate or accelerate with respect to the synchronously rotating air gap mmf, a relative 
motion begins. The equation describing the relative motion is known as the swing equation. 

Synchronous machine operation: 

• Consider a synchronous generator with electromagnetic torque Te running at synchronous speed cosm.. 

• During the normal operation, the mechanical torque Tm = Te. 

• A disturbance occur will result in accelerating/decelerating torque Ta=Tm-Te (Ta>0 if accelerating, Ta<0 if 
decelerating). 

• By the law of rotation - 

7 — = T =T —T 

Where J is the combined moment of inertia of prime mover and generator . 

• 6m is the angular displacement of rotor w.r.t. stationery reference frame on the stator 

• 6m = cosmt+5m, cosm is the constant angular velocity . 

• Swing equation in terms of inertial constant M 

• Relations between electrical power angle 5 and mechanical power angle 5m and electrical speed and 
mechanical speed 



where P is the pole number 
Swing equation in terms of electrical power angle 5 



2 d^-S 

- M = P - P 

p at- 



Con verting the swing equation into per unit system 



2Hd-8 _ _ 
^^dt^ ^n(fj:^j ^(pzO 



where M = — 
where H is the inertia constant . 

III. STEADY STATE STABILITY 

The ability of power system to remain its synchronism and returns to its original state when subjected to small 
disturbances. Such stability is not affected by any control efforts such as voltage regulators or governor. 

3.1 Analysis of steady-state stability by swing equation 

Starting from swing equation: - 



Hd-S 
nf^dt 



= p , ^ — P ' ^ = P — P sino 



Introduce a small disturbance A 5 
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• Derivation is from 5=50+A5 

• Simplify the nonlinear function of power angle 5 

• Analysis of steady-state stability by swing equation 

• Swing equation in terms of A5 

H d-AS 

— + R.. cosSi^AS = 0 
TT/o ai- 

3.2 Damping torque: 

• Phenomena: when there is a difference angular velocity between rotor and air gap field, an induction torque 
will be set up on rotor tending to minimize the difference of velocities. 

• Introduce a damping power by damping torque : 

dS 

• Introduce the damping power into swing equation. 

• Characteristic equation: 

H d^Ad dAd 

+ D^— + = 0 



■jt/o dt dt 
d'-AS dA3 ^ 

ai- at 



• Solution of the swing equation : 



d-A3 dA6 

, ^ + 2^o:i„— 1- oj-^AS = 0 

at- at 



• Roots of swing equation : 



Ad = ^"^^ e-^^^^sin (oj^t + 9} 



6 = 6n,-^ 



AS 



= e-<'^--sm(io^t-\-a') 



l-D- 





■O- 



Fig. 3.1 Diagram for Steady State Problem 



IV. TRANSIENT STABILITY 

The transient stability studies involve the determination of whether or not synchronism is maintained 
after the machine has been subjected to severe disturbance. This may be sudden application of load, loss of 
generation, loss of large load, or a fault on the system. In most disturbances, oscillations are of such magnitude 
that linearization is not permissible and the nonlinear swing equation must be solved. 



4.1 NUMERICAL SOLUTION OF SWING EQUATION 

The transient stability analysis requires the solution of a system of coupled non -linear differential 
equations. In general, no analytical solution of these equations exists. However, techniques are available to 
obtain approximate solution of such differential equations by numerical methods and one must therefore resort 
to numerical computation techniques commonly known as digital simulation. Some of the commonly used 
numerical techniques for the solution of the swing equation are: 

• Point by point method 

• Euler modified method 

• Runge-Kutta method 

In our analysis, we have used Euler modified method and Point-by Point Method. 
> The swing equation can be transformed into state variable form as : 
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dt 

dAS nff^ 



at 



p 
H " 



> We now apply modified Euler's method to the above equations as below 



dd 
'di 



Aw"^ = Atoi^ j_Li where Aoj^ij^^ = Aojf.^-\- 



dAoj 



at 



At 



dAto 



dt 



ill — -- 



d6 



^ ill wnere 0*^ ^j^^ = — 

at 



Ao}, At 



> Then the average value of the two derivatives is used to find the corrected values. 



= + 



0^ 

dt 



dt 



At 



dt 



aAoj 

0- H — 77- 

ai 



0 



At 



This is illustrated in the Simulink design. 
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Fig.4.1 Diagram for Transient State Stability design 



V. MULTI-MACHINE SYSTEM 

Multi-machine system can be written similar to one-machine system by the following assumptions: 
• Each synchronous machine is represented by a constant voltage E behind Xd (neglect saliency and flux 
change) . 

Input power remain constant . 

Using pre-fault bus voltages, all loads are in equivalent admittances to ground . 
Damping and asynchronous effects are ignored 
5mech = 5 

Machines belong to the same station swing together and are said to be coherent, coherent machines can 
equivalent to one machine . 



Solution to multi-machine system: 

• Solving initial power flow and determine initial bus voltage magnitude and phase angle 

S'i }Qi 



> Calculating load equivalent admittance 



Pj-JQ, 



> Nodal equations of the system 
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> Electrical and mechanical power output of machine at steady state prior to disturbances. 
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> 



Classical transient stability study is based on the application of the three-phase fault 



> Swing equation of multi-machine system 



= P 

77/ □ at- 



-z 

J-1 



Ei \\Ej\\Yij\cos[eij + 6j) = P^i - P 



EL 



> 
> 



Yij are the elements of the faulted reduced bus admittance matrix 
State variable model of swing equation 



dS. 

~dt 
dial 



= Ao}a , i = 1, Jf ,Ti 
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Fig. 5.1 Diagram for Dual-machine Stability 



VI. SIMULATION 

For our simulation we have taken a Power System Network of a electrical company shown in fig. 5.1 . 
The values of load data, voltage magnitude, generation schedule and the reactive power limits for the regulated 
buses are suitably taken. Bus 1 of the system is carefully considered with a specified value and is taken as slack 
bus . The data for basemva, accuracy, mailer, bus data, line data and gendata are suitably taken in matrix form . 
The whole process of simulation is carried out by using MATLAB-7.0. 
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Fig. 6.1 Dual-machine Stability for Fault cleared at 

0.4Sec . 
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Fig. 6. 2 Dual-machine Stability for Fault cleared at 

O.SSec . 
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Fig. 6. 3 Dual-machine Stability for Fault cleared at 0.697Sec . 



VII. RESULT ANALYSIS 

> Phase angle differences, after reaching a maximum of 521=123.90 and 

> 53 1=62.950 will decrease, and the machines swing together. 

> Hence, the system is found to be stable when fault is cleared in 0.4 second. 

> The swing curves shown in figure show that machine 2 phase angle increases without limit. 

> Thus, the system is unstable when fault is cleared in 0.5 second. 

> The simulation is repeated for a clearing time of 0.45 second which is found to be critically stable. 



VIII. CONCLUSIONS 

It can be seen that transient stability is greatly affected by the type and location of a fault so that a 
power system analyst must at the very outset of a stability study decide on these two factors. Thus we saw that a 
two-machine system can be equivalently reduced to a one machine system connected to infinite bus bar. In case 
of a large multi-machine system, to limit the computer memory and time requirements, the system is divided 
into a study subsystem and an external subsystem. The study subsystem is modelled in details whereas 
approximate modelling is carried out for the rest of the subsystem. The qualitative conclusions regarding system 
stability drawn from a two-machine or an equivalent one-machine infinite bus system can be easily extended to 
a multi-machine system. It can be seen that transient stability is greatly affected by the type and location of a 
fault so that a power system analyst must at the very outset of a stability study decide on these two factors. The 
solution of swing equation in Small Signal Stability Including Effects of Rotor Circuit Dynamics in synchronous 
machine will be the future work. 
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